Trajectories of mental health in children and adolescents during the COVID-19 pandemic: findings from the longitudinal COPSY study

Background Mental health and health-related quality of life (HRQoL) in children and adolescents deteriorated during the COVID-19 pandemic. The aim of this population-based longitudinal study was to explore whether distinct mental health trajectories in youths can be identified over the course of the pandemic. Methods Mental health problems (MHP), psychosomatic symptoms and HRQoL were assessed at five time points between May 2020 and October 2022 in 744 children and adolescents aged 7 to 20 years using established instruments. We used generalized mixture modeling to identify distinct mental health trajectories and fixed-effects regressions to analyse covariates of the identified profiles of change. Results We found five distinct linear latent trajectory classes each for externalising MHP and psychosomatic symptoms and four trajectory classes for internalising MHP. For HRQoL, a single-class solution that indicates a common development process proved to be optimal. The largest groups remained almost stable at a low internalising and externalising symptom level (64 to 74%) and consistently showed moderate psychosomatic symptoms (79%), while 2 to 18% showed improvements across the pandemic. About 10% of the youths had consistently high internalising problems, while externalising problems deteriorated in 18% of youths. Class membership was significantly associated with initial HRQoL, parental and child burden, personal resources, family climate and social support. Conclusions The mental health of most children and adolescents remained resilient throughout the pandemic. However, a sizeable number of youths had consistently poor or deteriorating mental health. Those children and adolescents need special attention in schools and mental health care. Supplementary Information The online version contains supplementary material available at 10.1186/s13034-024-00776-2.


Background
There is strong empirical evidence that the COVID-19 pandemic has posed substantial mental health challenges to children and adolescents.Systematic reviews and international surveys report considerable impairments of children's and adolescents' mental health during the pandemic [1][2][3][4][5][6].However, some single studies also show mixed results depending on the time of assessment, instruments, country and study design [6,7].Data suggests that both COVID infection dynamics and protection measures are positively associated with psychopathology [6,[8][9][10].Studies showing a deterioration of youth mental health mainly focus on the first 1.5 years of the pandemic [5,6].Most of those studies compare the pre-pandemic mental health with that during the pandemic or report two measurement points during the pandemic [7], but longitudinal studies are still scarce.
To our knowledge, the German nationwide longitudinal COPSY (COvid-19 and PSYchological Health) study is the first study exploring the mental health of children and adolescents across a time span of three pandemic years.In year one of the pandemic, general mental health problems, anxiety and depressive symptoms as well as the number of children with a low health-related quality of life (HRQoL) have been elevated compared to pre-pandemic data [11].In year two, mental health stabilised on an impaired level [12], while in year three an improvement became visible for HRQoL, general mental health problems and anxiety.Depressive symptoms even recovered to prepandemic levels [11].To our knowledge, there are hardly any comparable large-scale longitudinal published studies on the mental health of children and adolescents during the second and third year of the pandemic.
When reading the plentitude of studies published so far, it is surprising that most studies only report average results of their studied population.However, there is no empirical evidence on whether there are groups of children and adolescents who reacted similarly or differently to the pandemic challenges during the pandemic.Generalized mixture modeling (GMM) can be used to identify subgroups of participants with similar patterns of mental health trajectories.According to a recent review of studies, latent class analyses (LCA) and GMM are getting more common in child mental health research [13].However, we only found one paediatric study [14], which identified different mental health groups during the COVID-19 pandemic.Wang et al. [14] examined 2352 adolescents in China and found a resistance vs. a dysfunction group based on depressive and anxiety symptoms using latent growth mixture modeling of longitudinal data of three waves.Reviews of studies using GMM in adult data [15][16] summarise studies before and during the pandemic and report that approximately 66% of stress-exposed individuals respond resiliently, 13-21% show recovery responses, 9-12% show delayed mental distress, and about 11% report chronic distress in the face of general (and pandemic) stress.
In the present study, we aim to go a step beyond the previous statistics reported in most youth studies, which assume that the development in the study population follows a uniform model of change.We want to learn more about different mental health trajectories as well as potential predictors being associated with those trajectory groups.This allows a detailed characterisation of the development, which is essential to guide and target mental health prevention and intervention programs.We are particularly interested in exploring whether different group trajectories exist for internalising and externalising mental health problems and HRQoL because those are central constructs of positive and negative mental health, which have been heavily overlooked at the start of the pandemic and have been significantly impaired in youth during the pandemic.In the COPSY study those constructs were also chosen to allow comparisons to prepandemic studies.Further, previous studies showed that psychosomatic symptoms also increased during the pandemic but were widely overlooked and under-researched.Thus, we also added psychosomatic symptoms as a main construct to our analyses.
This study extends the previous analytical approach of the COPSY study, which was limited to the analysis of covariates of change.To understand potential predictors of trajectory groups of mental health, HRQoL and psychosomatic symptoms, we selected a variety of sociodemographic, personal, familial, and social factors.Previous publications found a gender effect in that girls showed lower HRQoL and more anxiety, depressive and psychosomatic symptoms during the pandemic than boys, while boys had more externalising mental health problems like hyperactivity and conduct behavior [11].Further, age effects were found in previous studies in that older females reported more anxiety [11].Those tendencies are in line with a recent review [6].Our COPSY study so far also found that youths with a good family climate, social support and personal resources were quite resilient [12,17] and other studies showed fewer negative effects or even benefits during the pandemic [18][19][20].Moreover, children may also have benefited from the pandemic, as they were able to spend more quality time with their parents during the lockdowns and may have avoided school and social problems [21,22].Further, the COPSY study identified important risk factors for child mental health and HRQoL.These include having a parent with a mental health problem or a parent that felt burdened by the pandemic, a migration background, low parental education and limited living space [11].Based on these previous research findings and their significance in the context of the study, gender, age, parental education, migration background, single parenting, parental pandemic burden, parental mental health problems, family conflicts, social support and personal resources were selected as predictors for the present analysis.
Overall, previous results of our COPSY study and the international literature led us to assume that there are at least two groups of different mental health trajectories.Thus, the objectives of this study are: 1. To explore whether distinct mental health, HRQoL and psychosomatic symptom trajectories in children and adolescents are identifiable over the course of the pandemic.We hypothesise that there is at least one large group of children and adolescents, who was resilient during the pandemic and showed either no impact or an initial decrease in mental health but recovered fast.Further, we hypothesise that there is at least one group characterised by decreasing mental health.2. To investigate predictors of mental health groups differing in their trajectories of improvement, stability or deterioration.We assume that age, gender, low parental education and migration background as well as single parents households, financial or other pandemic burden, previous mental health problems, family conflicts, low social support, and low personal resources may be associated with an increased risk of belonging to the identified low mental health group(s).

Study design and sample
The German population-based longitudinal COPSY study assessed children, adolescents and their parents at five time points: At the start of the pandemic, when there was a partial lockdown in Germany: 05-06/2020 (T1), during the first pandemic winter with a full nationwide lockdown: 12/2020-01/2021 (T2), after the summer in the second year, when infection rates were low and restrictions loosened: 09-10/2021 (T3), at the end of the second pandemic winter, when there were still regulations of private gatherings: 02/2022 (T4), and in autumn of year three, when only minimal restrictions remained: 09-10/2022 (T5).In each survey wave, between N = 1586 (T1) and N = 1701 (T5) families participated [11,23].
In this study, only the data of those n = 744 participants with complete data at all five measurement points were included who were aged 7 to 20 years and for whom parent-proxy report questionnaire data were available for all five measurement points.Respondents attending all five measurement points were older (1.3 years), were less likely to have a migration background and were more likely to have low parental education than respondents attending fewer measurements points.However, effect sizes were marginal.For single parenting, gender and occupation no statistically significant differences were observed.Regarding the mental health outcomes, a statistically significant difference was found for externalising mental health problems in three out of five measurements points, and for internalising mental health problems as well as for psychosomatic complaints in one out of five measurements points.Those participants attending all five measurements points displayed lower problems, however effect sizes again were only marginal.Families were invited via email to participate in the nationwide online survey which used quota sampling.This method helped to ensure that the sample reflected the sociodemographic characteristics of the German population.Detailed information on the study design and sample are published elsewhere [23,24].

Measures
The study included internationally established and validated instruments to assess the mental health of children and adolescents.HRQoL was measured via the parent-reported KIDSCREEN-10 Index [25].Means were calculated with higher values indicating higher HRQoL.Mental health problems were assessed using the parent-reported Strengths and Difficulties Questionnaire (SDQ) [26], which provides a total difficulties' score across 20 items and four subscales.Internalising (subscales emotional problems and peer problems) and externalising (subscales hyperactivity and conduct problems) scores were calculated, each ranging between 0 and 20, with higher scores indicating more severe problems.The frequencies of psychosomatic complaints like stomachaches, headaches, sleeping problems and irritability were measured via the parent-reported HBSC Symptom Checklist (HBSC-SCL) [27].The mean item score was calculated ranging from 1 to 5, with higher values indicating more psychosomatic complaints.All HRQoL and mental health measures (KIDSCREEN-10 Index, SDQ, HBSC-SCL) are internationally validated and have good psychometric properties, which have been described in detail previously [12].Familial resources were assessed using four items from the self-reported Family Climate Scale [28] and social support was measured using four items from the Social Support Scale [29].Higher mean values indicated more pronounced resources.
For sociodemographics, age and gender of the children and adolescents as well as their parents were assessed.The parents also reported their marital status, education, living space, and migration background.Children, adolescents, and their parents were further asked about their perceived burden of the pandemic using single items.

Statistical analyses
Mental health trajectories were explored using GMM to identify groups of children and adolescents, whose mental health changed in similar ways over the course of the pandemic.GMM is a person-centered approach, based on the assumption that the changes over the t = 5 measurement points in time can be adequately described if k distinct heterogeneous trajectory types underly [30][31][32][33][34][35][36][37].Each trajectory represents an underlying latent class.For all members of a specific latent class, the same linear development over time can be assumed.More specific, the GMM assumes that the observed data are generated by a mixture of multiple Gaussian normal distributions, each representing a distinct latent class (c j ; j = 1…k; unobserved subpopulations) [32,34].Different linear growth models are determined for each latent class.Within the GMM framework [37] the following assumptions were made: (1) Each person i belongs to one of k latent classes.
(2) Each of the k classes is characterised by the specific linear function y t|j = β 0j + β 0j • t s .The GMM approach uses a restricted maximum-likelihood approach to determine (1) the probability of belonging to one of the corresponding classes 1...k for each person i (p(cj |i)), and (2) the most likely regression parameter b 0j and b 1j for each class k [30,33].The estimation was done by determining the optimal values b 0j and b 1j under the assumption of j classes (1 < j < m).Based on the χ 2 -fit-value and the number of model parameters to be estimated, the information-theoretical Bayesian information criterion (BIC) was determined.The number of classes j, for which the BIC shows the optimal combination of data fit and parsimony (i.e.minimum BIC value), indicate the most plausible class solution [35].Additionally, the Lo-Mendell-Rubin adjusted Likelihood-Ratio (LMR-LR) test [31] indicated whether assuming j classes in contrast to j − 1 classes accounts for an incremental data fit: significant values indicated that the assumption of j classes is superior to the assumption of j − 1 classes.
Mplus was used to estimate the GMM [35].To avoid convergence in local maxima indicating invalid global solutions, which may have resulted in skewed, non-normal data distributions [30], Mplus used random starting values.Furthermore, the Expectation-Maximization-Algorithm was used in the iterative estimation process.
Bivariate associations of the latent classes with categorial variables of T1 were analysed by χ 2 -based techniques including the estimation of association strength by the contingency coefficient (CC).Class associations with continuous indicators were analysed by one-way analysis of variance and the effect sizes measure η (small: η = 0.1, medium: η = 0.25, strong: η = 0.37) [36].To quantify the respective change in each class between T1 and T5, a pairwise comparison of mean values was conducted using the effect size measure Cohen's d, with 0.2 being interpreted as a small, 0.5 as a medium and 0.8 as a strong effect.In order to be able to ensure the results of the variance-analytical testing of class differences in covariates with regard to distribution problems (e.g.floor effects), the tests were additionally examined using the non-parametric Kruskal-Wallis test.

Sociodemographics
The sociodemographic characteristics of the sample of n = 744 children and their parents are presented in Table 1.The age range of the children was between 7 and 20, with an initial mean of 12 years (SD = 3.3), the parents had a mean age of 44 years (SD = 7.1).There were slightly more boys reported upon than girls (51% vs. 49%).Most parents (59%) had a medium level of education, 14% had a migration background, 18% were single parents, 50% were full-time and 29% part-time employed.

Trajectories of mental health in children and adolescents
Mental health trajectories were investigated for internalising and externalising mental health problems, psychosomatic symptoms, and HRQoL.Using GMM, four to five classes were found for the former three outcomes-with the best data fit considering the number of estimated model parameters used.For HRQoL, the GMM showed an optimal data fit for the one-class solution, suggesting that a homogeneous change model valid for all can be assumed.Table 2 shows the class characterisations of internalising MHP, externalising MHP, and psychosomatic symptoms by associations to sociodemographics, risk and resource factors.

Internalising mental health problems
We found that four distinct mental health trajectory groups best described the data for internalising MHP (BIC: 17230.58,LMR-LR: 79.46, p =.002).The four groups are illustrated in Fig. 1 and described in The third (7%) and fourth class (9%) showed substantial improvements in internalising MHP (improvement groups) across the pandemic.

Externalising mental health problems
For externalising MHP, the best model fit of GMM analyses was found for a five trajectory groups solution (BIC: 16990.61,LMR-LR: 27.34, p =.184) as displayed in Fig. 2 and Table 3

Psychosomatic symptoms
With the GMM analyses a five psychosomatic symptom trajectory groups solution exhibited the most appropriate model fit (BIC: 3293.84,LMR-LR: 163.56, p =.007) (see Fig. 3; Table 3).The largest first class (impaired group, 79%) showed slightly worsening moderate symptoms (M T1 = 2.1 to M T5 = 2.3, Cohen's d = 0.28).This largest group can be characterised by a high HRQoL, the least child and parental burden, very good personal resources, the best family climate and high social support across the five classes.Interestingly there is a very small group (4%), which has a very similar but more stable trajectory (Cohen's d = 0.00) with the highest perceived child and parental burden and lowest personal resources of all those classes, yet a medium family climate.Parallel to those trajectories, two resilient groups (6 and 9%) showed stable trajectories (Cohen's d = 0.15/0.00)with no to very few psychosomatic symptoms.Finally, there is one class (improvement group, 2%), which initially had moderate psychosomatic symptoms, but show a steep decline of psychosomatic symptoms across the pandemic (Cohen's d = − 2.5).This group had a rather low initial HRQoL, was burdened the most initially and shows a medium child and parental burden, quite low initial personal resources and a bad family climate.Further findings of the covariate analyses indicated that there were few and rather weak associations between the class categorization of the three scales and most sociodemographic or COVID-related factors (see Table 2).Significant associations between the classes of each construct were found in the covariate analyses as displayed in Table 2 for parental burden (η = 0.193 to 0.290, p <.001), burden of the child (η = 0.224 to 0.298, p <.001), personal resources (η = 0.290 to 0.439, p <.001), family climate (η = 0.266 to 0.292, p <.001) and social support (η = 0.194 to 0.213, p <.001).
Class memberships proved to be highly correlated:

Discussion
This is one of the first studies using in-depth general mixture modeling of longitudinal mental health data over almost three years of the COVID-19 pandemic in children and adolescents.Therewith our COPSY study goes beyond reporting average statistics of mental health by investigating trajectories of children and adolescents.We identified four (internalising MHP) to five (externalising MHP and psychosomatic symptoms) distinct latent mental health trajectory classes.For HRQoL a single class solution was found to be appropriate.Thus, for HRQoL GMM analysis does not provide any further indications of informative differentiations of change processes that provide more detailed information on previously published COPSY findings.This was surprising to us and should be investigated in more detail in future studies.HRQoL can be defined as a global multidimensional construct which describes physical, social and psychological aspects of well-being and functioning [38], so it could be possible that it is too complex and global that different HRQoL trajectory groups could not be found.The present study found that almost three quarters of children and adolescents had consistently low internalising (64%) and externalising symptom levels (74%).Children belonging to this resilient group had good mental health and stayed mentally stable during the pandemic, confirming our first hypothesis that (at least) one group showed no pandemic impact.According to a prepandemic systematic review [13], a one large class solution with no or only minor symptoms is common when applying LCA/GMM to population-based mental health data.Our study confirms those findings for internalising and externalising MHP in times of health crises.A recent review on adult mental health studies using GMM [16] specified that, similar to a pre-pandemic meta-analysis [15], during the pandemic the pooled prevalence of the resilient groups was 66% in adults, which is between the 64% and 74% of resilient adolescents of the present study.
However, this pattern could not be found for psychosomatic symptoms.The majority (79%) of youths were allocated to the impaired group reporting moderate psychosomatic symptoms (once to four times per week), suggesting that psychosomatic reactions have been quite common since the start of the pandemic.This is in line with our previous COPSY publications [11].The most recent European HBSC study details that psychosomatic symptoms have been particularly frequent in 15-year-old girls with 50% of them reporting psychosomatic symptoms more than once a week during the pandemic [39].Another cross-national European study covering pre-pandemic data of children found that 46% had at least one frequent psychosomatic symptom during the past six months [40], which is similar to a review on adolescents and adults, which describes that 40 to 49% complained about at least one medically unexplained/ psychosomatic symptom [41].To summarise, findings of the present study indicate a strong increase in psychosomatic symptoms compared to the outlined large-scale studies and reviews.This might be explained by the fact that the COVID-19 infection itself led to psychosomatic complaints, that possible post-infection symptoms persisted or that children and adolescents psychosomatically expressed their worries or other emotions.However, evidence regarding these associations is limited and further studies are needed.
Moreover, the present study found three improvement groups for internalising and externalising MHP as well as psychosomatic symptoms.Thus, 2 to 18% of the initially moderately to severely impaired children and adolescents recovered over the pandemic.This also seems similar to the pooled prevalence of 13% allocated in recovery groups found in a review of adult studies during the pandemic using GMM [16].Our improvement groups found could be interpreted as meaning that the interests and needs of children and adolescents were heard and met after neglecting their mental health at the beginning of the pandemic, e.g., by adopting strict lockdown measures, so living conditions for these children may have improved.
Nevertheless, this study also revealed that a tenth of all children and adolescents (10%) have suffered from internalising MHP throughout the pandemic.This is very similar to a chronic adult mental health trajectory group (11%) found in pre-and pandemic adult studies using GMM according to reviews [15,16].In terms of externalising MHP, almost every fifth youth (18%) has had persistent or even increasing hyperactivity and conduct problems since the start of the pandemic.Those groups were called the poor mental health groups and the deterioration group, confirming the second part of our first hypothesis that there is at least one group that has suffered throughout the pandemic and is still burdened in the third year.This is in line with our previous research indicating that mental health problems increased during the pandemic and stabilised on an impaired level [11].These groups of children and adolescents need special attention in the context of prevention, diagnostics and mental health care.
The current study holds the advantage that it utilised GMM for an in-depth exploration of a longitudinal class characterisation of the sample.Many longitudinal pandemic paediatric studies using average statistics may oversimplify the complexity of intra-and interindividual variability [42].To our knowledge there are only few paediatric mental health studies, which used latent class modeling during pandemic times [43][44][45][46], most of which were cross-sectional studies.To our knowledge there is only one longitudinal adolescent study [47], which applied LCA in a similar manner as our COPSY study.Wang et al. [47] assessed mental health in a large Chinese adolescent's sample aged 16-25 years in 2020.They found two subgroups-called resistance vs. dysfunction-of adolescents based on their depressive and anxiety symptoms.In the whole sample they identified a subthreshold symptom increase in year 1, but when stratifying the sample in the two subgroups from February to June 2020, the larger resistance group (93%) had low depressive symptoms, which slightly increased but remained subthreshold, while the smaller dysfunction group (7%) had high depressive symptoms, which improved during the summer of 2020, but remained above the clinical cutoff score.The poor mental health groups found in our and the mentioned study may be related to risk factors like pre-existing mental or physical health problems, age, gender, a low socio-economic status, infection risk, lockdowns, low self-regulating abilities, low parental mental health or parenting quality, low family functioning, social support, isolation, loneliness, health-related worries or inconsistent routines (for a summary of risk factors, see [6]).
While in previous studies, some of the above-mentioned risk factors were found, this study surprisingly found that class membership was neither associated with any sociodemographic factors nor a COVID infection, which was contrary to our hypothesis 2-except for a small gender and age association for some groups of youths allocated in externalising MHP classes.
In line with hypothesis 2, we found that class membership was significantly associated with HRQoL, initial parental and child pandemic burden, personal resources, family climate and social support.Children and adolescents feeling burdened or with parents feeling burdened at the start of the pandemic felt worse during the pandemic than their peers.Youths with good personal resources, a good family climate and social support at the start of the pandemic were more frequently assigned to the resilient or improvement groups.This is in line with previous research on risk and protective factors, indicating that parental burden and psychopathology are associated with an increased risk of impaired mental health in children and adolescents during the pandemic, whereas personal, familial, and social resources support children and adolescents in staying mentally healthy [12,17].
The additional analyses regarding the associations between the class memberships of internalising, externalising and psychosomatic symptoms revealed high correlations between all classes of the constructs, indicating a strong content validity.The analyses found the largest association of the class members (80%) between the resilient internalising and resilient externalising classes and the largest association of class members of the psychosomatic impaired and both resilient internalising and externalising classes (86% and 75%).This could be due to the size of those classes which make up the majority of the youths analysed.The other associations are mixed and seem to be partly arbitrary, i.e. there were improved children in the same or the opposite classes of the other constructs with 31-61% overlap.
Findings of the present study, together with evidence from previous studies, suggest that more detailed analysis of mental health data is highly informative and worthwhile, as these methods help to expand our knowledge of interrelated factors in mental health and subgroups of mental health status at a given time point.For longitudinal data, trajectory modeling techniques like LCA and GMM or more individual-centered statistical approaches such as cluster analysis and sequence analysis produce important information on the evolution and progression of mental health subgroups over time, which is particularly needed in times of crisis like the COVID-19 pandemic.We would like to encourage researchers in the field to use these methods more frequently longitudinally to inform children's mental health care professionals and policymakers.
The strengths of this study are: (1) a large populationbased longitudinal sample of children and adolescents surveyed over a period of almost three years during the pandemic; (2) the use of internationally established and validated questionnaires; and (3) the application of advanced statistical methods.Limitations are: (1) that the results may not be generalizable to other populations or countries as the study was conducted in the context of the specific pandemic situation in Germany, which also differs from other countries worldwide in terms of wealth and aspects of healthcare (provision); (2) the SDQ and HBSC-SCL were used in this study, which, due to the item formulations and recorded content areas, exhibit floor effects in non-clinical, population-based samples.The characteristics of the identified classes must therefore be interpreted with regard to the rather demanding formulation ("minor problems" means minor problems understood in a clinical sense).It must be assumed that the dominance of the "resilient group" as the largest class is therefore also due to the low variability of the distribution of characteristics; (3) GMM is a data-driven, structure-identifying exploratory approach, for which replicative studies would be desirable for validation due to the lack of a priori hypothesis.Especially small classes could be statistical artifacts lacking reproducibility [48], thus we carefully selected the model using the fit indices described and conducted covariate analyses to substantiate their clinical meaningfulness.Although some of the identified classes are small, the identified class sizes are still above the critical limit of < 1% described in the literature.In addition to considering the class-specific fit indicators, the LMR-LR test was used to explicitly examine the need to consider the small classes in order to better and adequately describe the data information [49]; (4) there may be some content overlap between the investigated outcomes/scales (e.g. a few HRQoL items may be similar to SDQ internalising items), which may affect the discriminative interpretation of the results; (5) statistical power issues cannot be resolved in general terms, as the size of the classes to be identified is a decisive factor in this respect [49,50].For this reason, we chose an analysis approach that was as economic as possible, in that only those study participants who took part in all five measurement points were included in the data analysis.In addition, only linear trajectories were modeled, as estimating non-linear trajectories would have led to a disproportionate increase in the number of estimated parameters.The linearity assumption implies that only general stability and continuous improvement and deterioration during the pandemic could be identified; and (6) only those study participants who took part in all five measurement points were included in the data analysis.It cannot be ruled out that those excluded were not missing completely at random [51] and thus may bias the study results.

Conclusions
The mental health of children and adolescents has deteriorated during the COVID-19 pandemic, yet evidence on developmental trajectories based on longitudinal studies was limited.The population-based longitudinal COPSY study examined mental health trajectories and potential predictors associated with improvement, stability or deterioration.Most children and adolescents had good mental health and were resilient during the three pandemic years.However, one tenth to one fifth of the children had consistently high or deteriorating mental health problems, especially if they had low personal, familial and social resources, which call for special attention in schools and mental health care.Our findings provide mental health professionals and policymakers with evidence that can be used to support youth mental health in times of crises.

Fig. 1
Fig. 1 Trajectories of internalising mental health problems in children and adolescents during the COVID-19 pandemic

Fig. 3 Fig. 2
Fig. 3 Trajectories of psychosomatic symptoms in children and adolescents during the COVID-19 pandemic

Table 1
Sociodemographic characteristics of the COPSY sample M mean, SD standard deviation1Levels of education were operationalized based on the international 'Comparative Analysis of Social Mobility in Industrial Nations' (CASMIN) classification.This classification differentiates levels of education based on distinct combinations of academic and vocational qualifications.Based on these combinations, a categorization into parents with low (primary), medium (secondary), and high (tertiary) education was performed

Table 3
. The largest first class (resilient group, 74%) remained almost stable across the three pandemic years with very slight improvements at a low internalising symptom level (M T1 = 3.1 to M T5 = 2.7, Cohen's d = − 0.24).Participants belonging to this group showed the highest HRQoL, the lowest pandemic burden, high personal resources, the best family climate and best social support across all four classes.The second class (poor mental health group) with 10% of the participants started at a high symptom level at the beginning of the pandemic and remained at a high symptom level, improving only slightly (M T1 = 9.6 to M T5 = 8.8, Cohen's d = − 0.45).This group included participants with the lowest HRQoL, feeling most burdened by the pandemic and having the lowest personal resources and worst family climate compared to the other classes.
. Again, the largest first class (resilient group, 64%) remained almost stable with slight improvements at a low externalising symptom level (M T1 = 4.1 to M T5 = 2.9, Cohen's d = − 0.32) across the three pandemic years.Youths in that group tended to be more likely 12 than 10 years old, reported the highest HRQoL, the lowest child and parental burden, and the best personal resources, HRQoL and a low pandemic burden.Personal and familial resources were medium in that group.Further, two improvement groups (class 2 and 4, both 9%) showed a significant decrease (Cohen's d = − 0.76) in externalising symptoms across the five survey waves.And finally, the smallest group (deterioration group, 2%) starting at a high symptom level showed a substantial symptom increase across the pandemic (Cohen's d = 1.06).
family climate and social support compared to the other four classes.The second largest class (poor mental health group, 16%) showed increasing externalising symptoms from M T1 = 7.2 to M T5 = 8.7 (Cohen's d = 0.43), being rather 12 than 10 years old and having medium

Table 2
Associations between sociodemographics, risk and resource factors with the latent classes for internalising problems, externalising problems and psychosomatic symptoms 4%).Class 4 (impairment group) with regard to psychosomatic symptoms mainly includes members of internalising MHP class 3 (improvement group) (37.5%) and externalising MHP class 2 (improvement group) (40.6%).The contingency tables of the complete data distributions are attached in the Online Appendix A.